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Abstract 

In this paper, first we give a sequential linear-time algorithm for the longest 
path problem in meshes. This algorithm can be considered as an improve- 
ment of [l3|. Then based on this sequential algorithm, we present a constant- 
time parallel algorithm for the problem which can be run on every parallel 
machine. 
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1. Introduction 

The longest path problem, i.e. the problem of finding a simple path with 
the maximum number of vertices, is one of the most important problems in 
graph theory. The well-known NP-complete Hamiltonian path problem, i.e. 
deciding whether there is a simple path that visits each vertex of the graph 
exactly once, is a special case of the longest path problem and has many 
applications 0, 0] • 

Only few polynomial-time algorithms are known for the longest path 
problem for special classes of graphs. This problem for trees began with 
the work of Dijkstra around 1960, and was followed by other people 0, 0, 
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16l . Il8l . l22j. In the area of approximation algorithms it has been shown that 
the problem is not in APX, i.e. there is no polynomial-time approximation 
algorithm with constant factor for the problem unless P=NP [9j. Also, it has 
been shown that finding a path of length n — rf is not possible in polynomial- 
time unless P=NP [12| . For the backgrround and some known result about 
approximation algorithms, we refer the reader to [H, @, 0]. 

A grid graph is a graph in which vertices lie only on integer coordinates 
and edges connect vertices that are separated by a distance of once. A solid 
grid graph is a grid graph without holes. The rectangular grid graph R(n, m) 
is the subgraph of G°° (infinite grid graph) induced by V (m, n) = {v | 1 < 
v x < m , 1 < v y < n} , where v x and v y are respectively x and y coordinates of 
v (see Figure [T]). A mesh M(m, n) is a rectangular grid graph R(m, n). Grid 
graphs can be useful representation in many applications. Myers [19| suggests 
modeling city blocks in which street intersection are vertices and streets are 
edges. Luccio and Mugnia [17| suggest using a grid graph to represent a two- 
dimensional array type memory accessed by a read/write head moving up, 
down or across. The vertices correspond to the center of each cell and edges 
connect adjacent cells. Finding a path in the grid corresponds to accessing 
all the data. 

Itai et al. 



f lj have shown that the Hamiltonian path problem for general 



grid graphs, with or without specified endpoints, is NP-complete. The prob- 
lem for rectangular grid graphs, however, is in P requiring only linear-time. 
Later, Chen et al. |3| improved the algorithm of ll| and presented a parallel 
algorithm for the problem in mesh architecture. There is a polynomial-time 



algorithm for finding Hamiltonian cycle in solid grid graphs 15J. Also, the 
authors in 23|] presented sufficient conditions for a grid graph to be Hamilto- 
nian and proved that all finite grid graphs of positive width have Hamiltonian 
line graphs. 

Recently the Hamiltonian cycle (path) and longest path problem of a grid 



graph has received much attention. Salman et al. [21] introduced a family 
of grid graphs, i.e. alphabet grid graphs, and determined classes of alphabet 



grid graphs that contain Hamiltonian cycles. Islam et al. [10| showed that the 
Hamiltonian cycle problem in hexagonal grid graphs is NP-complete. Also, 
Gordon et al. [8| proved that all connected, locally connected triangular 
grid graphs are Hamiltonian, and gave a sufficient condition for a connected 
graph to be fully cycle extendable and also showed that the Hamiltonian 
cycle problem for triangular grid graphs is NP-complete. 

Moreover, Zhang and Liu 25] gave an approximation algorithm for the 
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Figure 1: The rectangular grid graph R(8, 7). 



longest path problem in gri d g raphs and their algorithm runs in quadratic 



time. Also the authors in [13[ has been studied the longest path problem 
for rectangular grid graphs and their algorithm is based on divide and con- 
quer technique and runs in linear time. Some results of the grid graphs are 



investigated in [14], |20 



In this paper, we present a sequential and a parallel algorithms for finding 
longest paths between two given vertices in rectangular grid graphs (meshes). 



Our algorithm has improved the previous algorithm [13[ by reducing the 
number of partition steps from 0[m + n) to only a constant. 

The organization of the paper as follow: In Section 2, we review some 
necessary definitions and results that we will need. A sequential algorithm 
for the longest path problem is given in Section 3. In Section 4, a parallel 
algorithm for the problem is introduced which is based on the mentioned 
sequential algorithm. Conclusions is given in Section 5. 

2. Preliminary results 

In this section, we give a few definitions and introduce the correspond- 
ing notations. We then gather some previously established results on the 
Hamiltonian and the longest path problems in grid graphs which have been 
presented in 0, [ill, 13 . 



The two-dimensional integer grid G°° is an infinite graph with vertex set of 
all the points of the Euclidean plane with integer coordinates. In this graph, 
there is an edge between any two vertices of unit distance. For a vertex v 
of this graph, let v x and v y denote x and y coordinates of its correspond- 
ing point (sometimes we use (v x ,v y ) instead of v). We color the vertices of 
the two-dimensional integer grid as black and white. A vertex v is colored 
white if v x + v y is even, and it is colored black otherwise. A grid graph G g 
is a finite vertex-induced subgraph of the two-dimensional integer grid. In 
a grid graph G g , each vertex has degree at most four. Clearly, there is no 
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Figure 2: A Hamiltonian cycle for the rectangular grid graph i?(5,4). 



edge between any two vertices of the same color. Therefore, G g is a bipartite 
graph. Note that any cycle or path in a bipartite graph alternates between 
black and white vertices. A rectangular grid graph R(m, n) (or R for short) 
is a grid graph whose vertex set is V(R) = {v | 1 < v x < m, 1 < v y < n}. 
In the figures we assume that (1,1) is the coordinates of the vertex in the 
upper left corner. The size of R(m, n) is defined to be mn. R(m, n) is called 
odd-sized if mn is odd, and it is called even-sized otherwise. In this paper 
without loss of generality, we assume m > n and all rectangular grid graphs 
considered here are oddxodd, evenxodd and evenxeven. R(m,n) is called 
a n-rectangle. 

The following lemma states a result about the Hamiltonicity of even-sized 
rectangular graphs. 

Lemma 2.1. 3/ R(m, n) has a Hamiltonian cycle if and only if it is even- 
sized and m, n > 1. 

Figure [2] shows a Hamiltonian cycle for an even-sized rectangular grid 
graph, found by Lemma 12.11 Every Hamiltonian cycle found by this lemma 
contains all the boundary edges on the three sides of the rectangular grid 
graph. This shows that for an even-sized rectangular graph R, we can always 
find a Hamiltonian cycle, such that it contains all the boundary edges, except 
of exactly one side of R which contains an even number of vertices. 

Two different vertices v and v' in R(m,n) are called color- compatible if 
either both v and v' are white and R(m, n) is odd-sized, or v and v' have 
different colors and R(m,n) is even-sized. Let (R(m,n), s,t) denote the 
rectangular grid graph R(m, n) with two specified distinct vertices s and t. 
Without loss of generality, we assume s x < t x . 

(R(m,n), s,t) is called Hamiltonian if there exists a Hamiltonian path be- 
tween s and t in R(m,n). An even-sized rectangular grid graph contains 
the same number of black and white vertices. Hence, the two end-vertices 
of any Hamiltonian path in the graph must have different colors. Similarly, 
in an odd-sized rectangular grid graph the number of white vertices is one 
more than the number of black vertices. Therefore, the two end-vertices of 
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any Hamiltonian path in such a graph must be white. Hence, the color- 
compatibility of s and t is a necessary condition for (R(m,n), s,t) to be 



Hamiltonian. Furthermore, Itai et al. 11] showed that if one of the following 



conditions hold, then (R(m,n), s,t) is not Hamiltonian: 

(Fl) R(m, n) is a 1-rectangle and either s or t is not a corner vertex (Figure 
Ha)) 

(F2) R(m,n) is a 2-rectangle and (s,t) is a nonboundary edge, i.e. (s,t) is 
an edge and it is not on the outer face (Figure 13(b)). 

(F3) R(m, n) is isomorphic to a 3-rectangle grid graph R'(m, n) such that s 
and t is mapped to s' and t' and all of the following three conditions 
hold: 

1. m is even, 

2. s' is black, t' is white, 

3. 4 = 2 and s' x < t' x (Figure He) ) or s' y ^ 2 and 4 < ^ - 1 (Figure 
Hd)). 



Also by [ll| for a rectangular graph i?(m, n) with two distinct vertices s and 
t, (R(m,n), s,t) is Hamiltonian if and only if s and t are color-compatible 
and R(m,n), s and £ do not satisfy any of conditions (Fl), (F2) and (F3). 
In the following we use P(R(m,n), s,t) to indicate the problem of finding a 
longest path between vertices s and t in a rectangular grid graph R(m,n), 
L(R(m,n), s,t) to show the length of longest paths between s and t and 
U(R(m,n), s,t) to indicate the upper bound on the length of longest paths 
between s and t. 

The authors in [13] showed that the longest path problem between any 
two given vertices s and t in rectangular grid graphs satisfies one of the 
following conditions: 

(CO) s and t are color-compatible and none of (Fl)- (F3) hold. 

(CI) Neither (Fl) nor (F2*) holds and either 

1. R(m, n) is even-sized and s and t are same-colored or 

2. R(m, n) is odd-sized and s and t are different-colored. 

(C2) 1. R(m, n) is odd-sized and s and t are black-colored and neither 
(Fl) nor (F2*) holds, or 
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Figure 3: Rectangular grid graph in which there is no Hamiltonian path between s 
and t. 



2. s and t are color-compatible and (F3) holds. 
Where (F2*) is denned as follows: 

(F2*) R(m,n) is a 2-rectangle and s x = t x or (s x = t x — 1 and s y ^ t y ). 

They also proved some upper bounds on the length of longest paths as fol- 
lowing: 

't x -s x + l, if(Fl), 

max{t x + s x , 2m - t x - s x + 2), if (F2*), 

mn, if (CO), 

mn-1, if (CI), 

mn-2, if(C2). 



U(R(m, n), s, t) 



Theorem 2.1. [IS] Let U(R(m,n), s,t) be the upper bound on the length 
of longest paths between s and t in R(m,n) and let L(R(m,n), s,t) be the 
length of longest paths between s and t. In a rectangular grid graph R(m, n), 
a longest path between any two vertices s and t can be found in linear time 
and its length (i.e., L(R(m,n), s,t)) is equal to U(R(m,n), s,t) . 



3. The sequential algorithm 

In this section, we present a sequential algorithm for finding a longest 
path between two vertices in rectangular grid graphs. This algorithm is the 
base of our parallel algorithm which is introduced in Section 4. First, we 
solve the problem for 1-rectangles and 2-rectangles. 

Lemma 3.1. fis}] Let P(R(m, n), s, t) be a longest path problem with n = 1 
or n = 2, then L(R(m, n), s, t) = U(R(m, n), s, t). 
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Figure 4: (a) Longest path between s and (ina 1-rectangle, (b) Longest path between s 
and t in a 2-rectangle, (c) and (d) A path with length 2m and 2m — 1 for a 2-rectangle, 
respectively. 



Proof. For a 1-rectangle obviously the lemma holds for the single possible 
path between s and t (see Figure 11(a)). For a 2-rectangle, if removing s 
and t splits the graph into two components, then the path going through 
all vertices of the larger component has the length equal to U(R(m,n),s,t) 
(see Figure 0(b)). Otherwise, let s be the vertex adjacent to s and t be the 
vertex adjacent to t such that s y ^ s y and t ^ t y . Then we make a path 
from s to s and a path from t to t as shown in Figure H(c), (d), and connect 
s to t by a path such that at most one vertex remains out of the path as 
depicted in this figure. □ 

From now on, we assume that m > n > 2, so one of conditions (CO), 
(CI) and (C2) should hold. Following the technique used in |3j we develop 
an algorithm for finding longest paths. 



Definition 3.1. [13( A separation of a rectangular grid graph R is a parti- 
tion of R into two disjoint rectangular grid graphs R\ and R 2 , i.e. V(R) = 
V(Rx) U V(R 2 ), and V{R X ) n V(R 2 ) = 0. 

Definition 3.2. llj Let v and v' be two distinct vertices in R. If v x < 2 
and v' x > m — 1, then v and u' are called antipodes. 

Definition 3.3. [3] Partitioning a rectangular grid graph R into five disjoint 
rectangular grid subgraphs R\ — R§ that is done by two horizontal and two 
vertical separations are called peeling operation, if the following two condi- 
tions hold: 

1. s,t G -R5 and s and t are antipodes. 
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2. Each of four rectangular grid subgraphs Ri — R4 is an even-sized rect- 
angular grid graph whose boundary sizes are both greater than one, or 
is empty. 

Generally the two vertical separation of a peeling are done before the two 
horizontal separation. However, for an odd x odd or odd x even rectangular 
grid graph with s x = t x , this order is reversed in order to guarantee that the 
boundary sizes of R3 and R4 are greater than one. Figure [5] shows a peeling 
on R(15, 11) where s is (6, 5) and t is (8, 9). 

The following lemma can be obtained directly from Definition 3.3. 
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Figure 5: A peeling on i?(15, 11) 



Lemma 3.2. Jj/ Let i?5(n 5 , m 5 ) be the resulting rectangular grid subgraph of 
a peeling on R(n,m), where s,t G V(R 5 ). Then 

1. s,t remain the same color in R 5 as in R; and 

2. R§ has the same parity as R, that is, 7715 mod 2 = m mod 2, and n§ 
mod 2 = n mod 2. 

Definition 3.4. A peeling operation on R is called proper if \R\\ + I-R2I + 
I-R3I + \R±\ + U(R 5 ,s,t) = U(R(n,m), s,t), where \Ri\ denotes the number 
of vertices of 

Lemma 3.3. For the longest path problem P(R(n,m),s,t), any peeling on 
R(m, n) is proper if: 

1. The condition (CO) holds and m mod 2 = n mod 2 {i.e. R(m,n) is 
even x even or odd x odd), or 
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2. One of the conditions (CI) and (C2) hold and R(m, n) is even x odd 
or odd x odd. 



Proof. The lemma has been proved for the case that (CO) holds (see |3j). So, 



we consider conditions (CI) and (C2). From Lemma [3T2l we know that s and 
t are still color-compatible, and we are going to prove that P(R 5 (m 5 , n 5 ), s, t) 
is not in cases Fl and F2*. 

By Lemma EJT2"} when R(m, n) is an oddx odd rectangular grid graph, R 5 (m^, n$) 
is also an odd x odd rectangular grid graph, s and t have the same color as 
in R, and hence R^{m^,n^) is not a 2-rectangle. If R$ is a 1-rectangle, then 
s y = ty or s x = t x and then we have the two following cases: 

Casel. (CI) holds and both s and t are different color. In this case, one 
of s x and t x (s y and t y ) is even and the other is odd. Considering that s and 
t are antipodes and R 5 is oddx odd, one of s and t must be at the corner and 
exactly one of the vertex goes out of the path. 

Case2. (C2) holds and both s and t are black color. In this case, all s x , 
s y , t x and t y are even. Hence, vertices s and t are before corner vertices and 
exactly two vertices go out of the path. 

In the similar way, when R(m, n) is an even x odd rectangular grid graph 
((CI) holds), i?5(m5,n 5 ) is also an even x odd rectangular grid graph, and 
hence R^{m^,n 5 ) is not a 2-rectangle. If R$ is a 1-rectangle, then s y = t y . 
In this case, s x and t x are both odd or even. Hence, s or t are at the corner 
and exactly one vertex goes out of the path. 

Therefore by Theorem 12. 1\ U(R 5 (m 5 ,n 5 ),s,t) = U(R(m,n), s,t) and any 
peeling of R(m, n) is always proper. □ 

Nevertheless, a peeling operation in an even x even rectangular grid graph 
R(m,n) may not be proper, and U(R 5 (m 5 ,n 5 ), s,t) ^ U(R(m,n), s,t), see 
Figure E] where the dotted-lines represent a peeling operation. In the two 
following peeling operation is not proper: 

(Fl') s is black, s x is even (or odd), t y = s y + 1 and s x ^ t x 
(F2') s is white, s x is even (or odd), t y = s y — 1 and s x ^ t x ; 

Lemma 3.4. For the longest path problem P(R(n,m), s,t) , where R(m,n) 
is an even x even rectangular grid graph, a peeling operation on R(m, n) is 
proper if and only if P(R(n, m), s, t) is not cases in (Fl ) and (F2 ). 

When a peeling operation is not proper it can be made proper by adjustment 
the peeling boundaries. In that case, if R lf R 2 , R3 and i? 4 are empty, then 
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Figure 6: Rectangular grid graph in which a peeling operation is not proper. 
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Figure 7: 



('1) 



# 5 is 2-rectangle that is in case (F2*). Therefore, without loss of generality, 
we assume R\, R 2 , R3 or #4 is not empty. If rectangular grid subgraphs #3 
and R4 are empty, then we move one column (or two columns when R\ or R2 
is a 2-rectangle) from Ri or R 2 to #5 such that #1 or R 2 is still even-sized 
rectangular grid graphs; see Figure [7^ a). If #1, i? 2 , # 3 and # 4 (or # 3 and 
R4) is not empty, then we move one row (or two rows when R 3 or #4 is 
2-rectangle) from i? 3 or #4 to #5 (Figure ITJb)), or move the bottom row to 
i?4 (Figure [7(c)) or move the upper row to f? 3 (Figure UJd)), such that #3 or 
i?4 is still even-sized rectangular grid graphs. 

After a peeling operation on #(771,71), we construct longest paths in 
#5(7715,715). Consider the following cases for #5(7715,715): 

(a) 7715, 715 < 3. 

(b) 7715, 715 are even, and either 7715 > 4 or 715 > 4; 

(c) 77i 5 , n 5 are odd, and either 7715 > 5 or ti 5 > 5; 

(d) 7715 is even and ti 5 is odd, and either 7715 > 4 or ti 5 > 5. 

For case (a), we showed that when 71 = 1, 2 the problem can be solved easily. 
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(c) 

Figure 8: For n = m = 3, (a) s and t are white, then there is Hamiltonian path, 

(b) s and t have different colors, then there is a path with U(R, s, t) = mn — 1 and 

(c) s and t are black, then there is a path with U(R, s, t) = mn — 2. 

For m, n = 3 the longest paths of all the possible problems are depicted in 
Figure M (the isomorphic cases are omitted). 

For cases (b), (c) and (d) we use the definition of trisecting. 

Definition 3.5. [i[ Two separations of R$ that partition it into three rect- 
angular grid subgraphs Rl, R\ and R™ is called trisecting, if 

(z). i?5 and are a 2-rectangle, and 

(ii). s G V(R S 5 ) and t G V{R\). 

A trisecting can be done by two ways horizontally and vertically. If 
< 4 or m^,n^ > 4, then trisecting is done horizontally, if 715 < 4, then 
trisecting is done vertically. 

Definition 3.6. A corner vertex on the boundary of R s b (resp., R\) facing 
R™ is called a junction vertex of R^ (resp., R\) if either 

(i) The condition (CO) holds and it has different color from s and t, or 

(ii) One of the conditions (CI) or (C2) hold and U(R S 5 , s,p)+U(R^, m, m')+ 
U(R\, q, t) = U(R(m, n), s, t). Where p is one of the corner vertices of 
i?5, q is one of the corner vertices of R\, and m and m are two of the 
corner vertices of R™ facing R s b and R\, respectively. 
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Figure 9: A trisecting on i?(6, 6). 



In Figure pi and 77 2 , ?i and 52, wi, 777,2, 7773 and m 4 are junction vertices 
in i?5 and R™, respectively. Existence of junction vertices has been 
proved for condition (CO) in [3|, in this paper we only consider conditions 
(CI) and (C2). 

Lemma 3.5. Performing a trisecting on R^,, where m,n > 3, assuming 
condition (CI) or (C2) holds, if n 5 = A, and s and t facing the common 
border of R^ and R\, then there is no junction vertex for R^ and R\, otherwise 
i?5 and R\ have at least one junction vertex. 

Proof. Consider Figure ITOT a) and (b), where n 5 = 4 and two vertices s 
and t facing the common border R s 5 and R\. In this case, the only two 
vertices p 2 and q\ may be junction vertices. By Theorem I2.1[ there exists a 
Hamiltonian path from s to p 2 and from qi to t in R^ and R\, respectively, 
and U(Rl, s,p 2 ) + U(R t 5 ,q ll t) 7^ U{R(m,n), s,t). Hence, neither R^ nor R\ 
has a junction vertex at all. Now for other cases, we show that R s b and R\ 
have at least one junction vertex. 

In case (b), s and t have the same color (white or black), and two corner 
vertices on the boundary of R s b (resp., R\) facing R™ are different color and 
also R™ is a k— rectangle that k is empty or k > 2 and even. We consider 
the following three cases for s and t: 

Case 1. Both s and t are the corner vertices on the boundary of R^ and 
R\ facing R™; see Figure [TDT c). By Theorem 12 .![ there exists a Hamiltonian 
path from s to p 2 and from qi to t in R^ and R\, respectively, and a path 
from 777,3 to 777-2 which does not pass through a vertex in R™. Therefore, 
U(Rl, s,p2) + U(R£ l , 777 3 , m 2 ) + U(R\, qi,t) = U(R(m, 77), s, t) and hence both 
i?5 and R\ have a unique junction vertex. 

Case 2. s is the corner vertex on the boundary of R s b facing R™; see Figure 
ITOT d). By Theorem I2.1[ there exists a Hamiltonian path from s to p 2 , from 
qi to t and a path from 7773 to 7772 which does not pass through a vertex, or a 
Hamiltonian path from s to p 2 , form 777 2 to 7774 and a path from q 2 to t which 
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Figure 10: A trisecting on i?(6,6), i?(6,4). 



does not pass through a vertex. Therefore, U{R%, s,p 2 ) + U(R™, m 3 , 777,2) + 
U{Rl, qi ,t) = U(R(m 1 n) 1 s,t) a ndU(Ri 1 s,p 2 )+U(R^,m 2l m A )+U(R t 5 ,q 2 ,t) = 
U(R(m, n), s, £) and hence i?g has a unique junction vertex and R\ have two 
junction vertices (the same argument is also applied to £). In this case, where 
n 5 = 4, both i?5 and R\ have a unique junction vertex; see Figure HUlfe). 

Case 3. s and £ are not the corner vertices on the boundary of R s 5 and 
R\ facing ft!™; see Figure [TDT f). By Theorem 12.11 there exists a Hamilto- 
nian path from s to p\, mi to m 3 and a path from q± to £ which does not 
pass through a vertex, or a Hamiltonian path from s to pi, from g 2 to £ and 
a path from mi to 7774 which does not pass through a vertex, or a Hamil- 
tonian path form 7772 to 7774, from q 2 to £ and a path from s to ^2 which 
does not pass through a vertex. Therefore, U{Ri ) , s,pi) + U(R$, mi, m 3 ) + 
U(Rlqi,t) = U(R{m,n),s,t), U(R s 5 ,s,pi) + U(R™,mi,m 4 ) + U{R\,q 2 ,t) = 
U(R(m,n),s,t) and U{Rl,s,p 2 )+U{R£,m 2 ,m 4 )+U{Rl,q 2 ,t) = U{R{m,n),s, 
and hence both R^ and Rl have two junction vertices. 

In case (c), s and £ are black or different color, and two corner vertices 
on the boundary of R^ (resp., RQ facing R™ are black and also R™ is a 
A;— rectangle that k > 1 and odd. There are three cases for s and £: 

Case 1. Both s and £ are the corner vertices on the boundary of R^ and 
R\ facing f?™; see Figure UlTa). Then s and £ are black. By Theorem 12.11 
there exists a path from s to p 2 , from 51 to £ which does not pass through a 
vertex in R^ and R\, respectively, and a Hamiltonian path from 7773 to 777,2 in 
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Kg. Therefore, U{R%, s, p 2 ) + U(Rf, m 3 , m 2 ) + U{R\, q u t) = U(R(m,n),s,t) 
and hence both R^ and R\ have a unique junction vertex. 

Case 2. s is the corner vertex on the boundary of R s b facing R™, then s is 
black and t is black or white; see Figure [TTT b). By Theorem 12 .![ there exists a 
path from s to p 2 and from qi (or q 2 ) to t, where t is black, which does not pass 
through a vertex, and a Hamiltonian path from and m 2 to m 4 (or from 777.3 to 
7772), or a Hamiltonian path from gi (or q 2 ) to £, where t is white and 777-2 to m 4 
(or from 777,3 to 7772) and a path from s to p 2 which does not pass through a ver- 
tex . Therefore, U(R s 5 ,s,p 2 ) + U(Rg,m 3 ,m 2 ) + U(Ri,q 1 ,t) = U(R(m,n), s,t) 
and U(Rl, s,p 2 ) + U(R™, m 2 , 7774) + U(R\, q 2 ,t) = U(R(m, 77), s, £) and hence 
i?5 has a unique junction vertex and R\ have two junction vertices (the same 
argument is also applied to t). In this case, where n 5 = 5, both R^ and a 
unique junction vertex; see Figure [TTT c). 

Case 3. s and t are not the corner vertices on the boundary of i?g and 
R\ facing R™\ see Figure ITlT d). By Theorem 12. 1[ there exists a Hamiltonian 
path from s to p, m to 777 and g to t, where s (or t) is white, and a path form 
s to p and 9 to t which does not pass through a vertex where s (or t) is black, 
jo is pi or 772, (7 is q\ or g 2j 771 is mi or 7772 and m is 7773 or 7774. Therefore, 
U(R s 5 ,s,p) + U(Rf,m,m) + U(Rl,q,t) = U(R(m,n), s,t) and hence both 
i?5 and have two junction vertices. 

In case (d), if 775 > 3, the trisecting is performed horizontally, and the claim 
is proved by applying the same argument for case (b); see Figure UTT e). If 
77 5 = 3, the trisecting is performed vertically and also two corner vertices on 
the boundary of R^ facing R™ are black. Therefore, the claim is proved by 
applying the same argument for case (c). □ 

After trisecting, we construct a longest path in R^, R™ and R\ between 
s and p, 777 and 777 and q and t, respectively. In the case that none of R^ and 
R\ have junction vertices (when 775 = 4 and both s and t facing the common 
border R s b and i2|), we construct a longest path in R^ (resp., R\) between 
s (resp., t) and a none-corner vertex of the boundary facing R\ (resp., R^); 
see Figure H2J At the end, the longest paths in R 5 are combined through the 
junction vertices; see Figures [TU1 and [TT1 

Then we construct Hamiltonian cycles in rectangular grid subgraphs R\ 
to i?4, by Lemma I2.1( see Figure [13j Then combine all Hamiltonian cycles 
to a single Hamiltonian cycle. 

Two non-incident edges ei and e 2 are parallel, if each end vertex of ei 
is adjacent to some end vertex of e 2 . Using two parallel edges e\ and e 2 of 
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(d) (e) 

Figure 11: A trisecting on R(7,7), R(7,5) and R(6,5). 

two Hamiltonian cycles (or a Hamiltonian cycle and a longest path), such as 
two darkened edges of Figure [T4l(a). we can combine them as illustrated in 
Figure [TU^b) and obtain a large Hamiltonian cycle. 

Combining the resulted Hamiltonian cycle with the longest path of R§ is 
done as in Figure [15] 



nrun t 




Figure 12: the Longest path in i?(6,4). 

Considering all of the above, we get the algorithm of finding a longest 
path in rectangular grid graphs, as shown in Algorithm 13.11 

Consider the pseudo-code of our algorithm in Algorithm 13.11 The step 1 
dose only a constant number of partitioning, during the peeling operation, 
which is done in constant time. The step 2 trisects R$ which requires also a 
constant number of partitioning. Then finds a longest path in i?5 by merging 
paths of the partitions which can be done in linear time. The step 3 finds 
Hamiltonian cycles of Ri to R4 which is done in linear time. The step 4 
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Figure 13: Hamiltonian cycles in R\ to R4. 
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(a) (b) 
Figure 14: Combining two Hamiltonian cycles. 



which combines the Hamiltonian cycles and the longest path requires only 
constant time. Therefore, in total our sequential algorithm has linear-time 
complexity. 



4. The parallel algorithm 

In this section, we present a parallel algorithm for the longest path prob- 
lem. This algorithm is based on the sequential algorithm presented in the 
previous sections. Our parallel algorithm runs on every parallel machine, 
we do not need any inter-processor connection in our algorithm. We assume 
there are nm processors and they work in SIMD mode. For simplicity, we use 
a two-dimensional indexing scheme. Each vertex v of the given rectangular 
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Figure 15: The longest path between s and t. 



Algorithm 3.1 The longest path algorithm 
procedure LongestPath (R(m, n), s, t) 

Step 1. By a peeling operation, R(n, m) partitions into five disjoint rectangular 

grid subgraphs R\ to R§, such that s,t € R§ 

Step 2. Finding longest path between s and t in R$. 

Step 3. Construct Hamiltonian cycles in rectangular grid subgraphs i?i to R^. 

Step 4. Construct a longest path between s and t by combine all Hamiltonian 
cycles and a longest path. 



grid graph R(m,n) is mapped to processor (v x ,v y ). Each processor knows 
its index, coordinates s and t, and m and n. 

The peeling phase is parallelized easily, every processor calculates the 
following four variables, in parallel 
s x — 2; s x mod 2 = 
s T — 1; otherwise 



t x + 1; t x mod 2 = m mod 2 
t x + 2; otherwise 

min(s y ; t y ) — 2; mm(sj / ; t^) mod 2 = 
min(s y ; t y ) — 1; otherwise 

max(s y ] t y ) + 1; maa;(sj / ; t y ) mod 2 = n mod 2 
max{s y ] t y ) + 2; otherwise 

Where variables r 1; r 2 , r 3 and r 4 correspond to the right-most column 
number of R\ , the left-most column number of i?2 , the bottom row number of 



r 3 



r 4 
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i?3, and the top row number of R4, respectively. Then a processor can identify 
its subrectangular by comparing its coordinates with these four variables. 
In case (Fl ) and (F2 ), the boundary adjustment can be done by simply 
decrementing R\, R2, R3 or R4 or incrementing R3 or R±. 

The trisecting phase is also parallelized in a similar manner. In the fol- 
lowing we describe how we parallelized the horizontal trisecting, in two cases 
when R(m,n) is evenxodd (or oddxodd) and when it is evenxeven. In case 
R(m,n) is evenxodd or oddxodd, every processor simultaneously calculate 
the following two variables: 



Where variables I and r correspond to the bottom row number of R§ 
(resp. and the top row number of R\ (resp. Rl), respectively. 

In case R(m,n) is evenxeven, every processor simultaneously calculate 
the following two variables: 



A similar method can be used to parallelize the vertically trisecting. 
After peeling and trisecting, all processors in the same subrectangles simulta- 
neously construct either a longest path, Hamiltonian path or cycle according 
to the pattern associated with the subrectangle. For constructing a Hamil- 
tonian path in a rectangular grid graph, we use the constant-time algorithm 
of (if. For constructing a Hamiltonian cycle in an even-sized rectangle, we 
use the constant-time algorithm of [4] in which every processor computes its 
successor in the cycle. This algorithm is given in Algorithm 14. 1( see Figure 



For constructing a longest path, parallel algorithms can be easily devel- 
oped for each different pattern shown in Figure [ToT b). (c). As two examples, 
for constructing a longest path between vertices (2,1) and (m — l,n) in an 
oddxodd rectangular grid graph R(m, n), and vertices (m— 1, 1) and (n—1, 1) 
in an evenxeven rectangular grid graph R(m, n). We have developed the sim- 




min(s y ; t y ) min(s y ; t y ) mod 2 = 

min(s y ; t y ) + 1 min(s y ; t y ) mod 2^0 

max(s y ; t y ) max(s y ; t y ) mod 2 = 

max(s y ; t y ) — 1 max(s y ; t y ) mod 2^0 




min(s y ; t y ) min(s y ; t y ) mod 2 = 

min(s y ] t y ) + 1 m.in(s y ; t y ) mod 2^0 

max(s y ; t y ) max(s y ; t y ) mod 2^0 

max(s y ; t y ) — 1 max(s y ; t y ) mod 2 = 



da) . 
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pie algorithms Algorithm 14.21 and I4.3[ respectively. The algorithms for other 
patterns can be derived in the similar way. 

Then combining phase is parallelized as follows. The two processors at 
the two endpoints of a corner edge in a Hamiltonian cycle c\ check whether 
a neighboring Hamiltonian cycle C2 exists or not. If c 2 exists, then their 
successors are changed to the adjacent processors in c 2 . Similarly, the two 
processors at the endpoints of a corner edge in the longest path P in R 5 also 
check the existence of the adjacent edge in the Hamiltonian cycle C, and 
change their successors. Thus, the combining phase can be parallelized in 
constant steps without inter-processor communication. 




(a) (b) (c) 

Figure 16: (a) A Hamiltonian cycle in i?(4, 7), (b) and (c) two patterns of longest path in 
i? 5 (5,5) and R 5 (8,6). 



Algorithm 4.1 The Hamiltonian cycle parallel algorithm for an even-sized 
rectangular grid graphs 

procedure LongestPathi?(m, n) 

1: for each processor (x,y) in R(m,n) do in parallel 

2: if y = 1, then successor (x,y) ^— (x + l,y) 

3: elseif (y = 2, x is odd and x ^ 1) or (y = n and x even), then successor 

(x,y) < — (x - l,y) 
4: slesif x is even and y < n, then successor (x, y) < — (x, y + 1) 
5: else x is odd and y < n, then successor (x,y) < — (x,y — 1) 



5. Conclusion and future work 



We presented a linear-time sequential algorithm for finding a longest path 
in a rectangular grid graph between any two given vertices. Since the longest 
path problem is NP-hard in general grid graphs ll( , it remains open if the 
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Algorithm 4.2 The longest path parallel algorithm for oddxodd rectangular 
grid graphs 

procedure LongestPath (R 5 (m 5 , n 5 ), s, t) 

1: for each processor (x,y) in .R5 (7775, 775) do in parallel 

2: if x = 1 and y = 1 or x = m and y = n, then successor (x, y) < — nill 

3: elseif y is odd and x < m, then successor (x,y) < — (x + l,y) 

4: slesif y is odd and x = m, then successor (x, y) < — (x, y + 1) 

5: elseif y is even and x > 1, then successor (x,y) < — (x — l,y) 

6: else y is even and x = 1, then successor (x,y) ^ — (x,y + 1) 



Algorithm 4.3 The longest path parallel algorithm for even x even rectan- 
gular grid graphs 

procedure LongestPath (i? 5 (m 5 , 775), s, t) 

1: for each processor (x,y) in ^5(7775,775) do in parallel 
2: if x = m and y = 1, then successor (x,y) ^ — nill 

3: elseif (y is odd and x = m), (y is even and x = 1) or (y = n and x is even) 

then successor (x, y) < — (x, y — 1) 
4: slesif (y is odd and x < m), {y = n — 1 and x is even), (y = n and x is odd) 

then successor (x, y) < — (x + 1, y) 
5: elseif y is even and x > 1 then successor (x,y) < — (x — l,y) 
6: elseif y = n — 1 and x is odd then successor (x, y) i — (x, y + 1) 
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problem is polynomially solvable in solid grid graphs. Based on the sequential 
algorithm a constant-time parallel algorithm is introduced for the problem, 
which can be run on every parallel machine. 
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